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Abstract. We present the results of a detailed analysis of the XMM-Newton observation of the galaxy cluster 
^ ' Abell 3921. The X-ray morphology of the cluster is elliptical, with the centroid offset from the brightest cluster 

galaxy by 17", and with a pronounced extension toward the NW. Subtraction of a 2D /3-model fit to the main 
cluster emission reveals a large scale, irregular residual structure in the direction of the extension, containing both 
*0 ' diffuse emission from the intra cluster medium, and extended emission from the second and third-brightest cluster 

galaxies (BG2 and BG3) . The greatest concentration of galaxies in the subcluster lies at the extreme northern edge 
of the residual. The cluster exhibits a remarkable temperature structure, in particular a bar of significantly hotter 
gas, oriented SE-NW and stretching from the centre of the cluster towards BG2 and BG3. Our detailed study of 
the morphological and thermal structure points to an off-axis merger between a main cluster and a less massive 
«4 ' galaxy cluster infalling from the SE. From comparison of the temperature map with numerical simulations, and 

with independent calculations based on simple physical assumptions, we conclude that the merging event is ~ 0.5 
Gyr old. The cluster is thus perhaps the best X-ray observed candidate so far of an intermediate mass ratio, 
S_i ' moderate impact parameter merger. 

Key words, galaxies: clusters: general - galaxies: clusters: individual: Abell 3921 - galaxies: clusters: intergalactic 
medium, mergers - X-rays:general - X-ray:galaxies:clusters cosmology: large-scale structure of Universe 



^ ! 1. Introduction tion of the ICM (e.g., Roettiger, Burns & Loken 

. . . Ricker 1998). Mergers have also been found to modify the 

In the current scenario for the formation of structure m , , , ,. ,& Jn . h two . 

member galaxy distribution (bchmdler & Bohrmger 1993). 

the Universe, clusters of galaxies form hierarchically by a mi ■ , /■ ,i T nii • n i j- ±. 

' ? i lhus, substructures (m the ICM or m the galaxy distnbu- 

succession of mergers with objects already formed. During ,. . , , r ■-, j r n v.- a 

° J fi , t . . i tion) m clusters are a tossil record ot trie merger history, 

a merger event more than 10 ergs are dissipated m the ... . . , . „ 

. . . ,. ,.„.., ,. , , , ,. , X-ray observations ot galaxy clusters lorman & 
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r . ... i i,.,. ,. in Jones 1982 Mohr, Fabncant & Gellcr 1993) suggest that 

lorm initially between the colliding objects and then prop- , , . , , , , - r , ^ . ,„ 

. .. . , .. ° , . such objects at low redshift are dynamically young (Buotc 

agate m the impact direction. Ihese merger events have „ m .. 1 AA a. .... , . r . , , 

„ A , , . . , , ,. & Isai 1996.), displaying a large degree ot density and tem- 

a strong enect on the cluster, triggering star formation , , , TT » V. • . u aaai » , , ., , 

.. . n, u . , , . ,° ... ° . , . nun . . perature structure (e.g., Henry & Briel 1996; Markevitch 

(EvrardHaSp and gas stripping (Fujita et a lffl mthe & vikhlinin 2001; V MazzQtt Ed & Markevitch 2003; 
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galaxies, and leading to strong variations oi the physical _. T , .. .Viu,. „. ., , , , 

n i i ■. i ,i Neumann et al. 2003). Since the temperature substruc- 

properties (such as the temperature, the density and the . „ ; ,. ,. , ., . , 

\ „ , , T , ture survives tor a longer time than density substructure 

entropy) of the ICM. , ± . . , c , . ° . x I 

. . (typically 4 to 6 times longer), the temperature is a strong 

rrom the theoretical point oi view, numerical N-body . ,. r , , . , , ., , n . , 

, . , , , , , indicator ot the cluster history and its present dynamical 

simulations have shown that mergers produce substruc- s ^- a ^- e 

ture observable in the density and temperature distribu- „..,,,. „ , . . , 
Statistical studies oi cluster morphology can thus pro- 
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of structure formation to statistical studies of nearby clus- 
ters has been until now limited by the knowledge of the 
effect of cluster formation and evolution on the observa- 
tional properties. Despite the large amount of data avail- 
able, the dynamical evolution of the 1CM and the relation 
between the galaxies and gas during a merger event is still 
poorly understood (see Buote l2*0"0"2l for a review). A better 
understanding of the physics of merger events, in particu- 
lar the relaxation time of substructures fNakamura ll995|) 
is required. 

Combined X-ray and optical studies (e.g. Roettiger 
et al. [H)M Durret et al. ITMHl Ferrari et al. have 
the advantage of combining information on the collisional 
component (the gas), and the nominally non-collisional - 
and thus more likely to trace the dark matter - component 
(the galaxies). Such studies are thus well suited to inves- 
tigate in deeper detail the dynamical process of merging 
clusters. 

In this paper we present the results obtained from 
the XMM-Newton observation of the merging cluster 
A3921. This is the second target from our XMM-Newton 
Guaranteed Time (GT) program, which was established 
as a systematic study of merging clusters of galaxies, with 
the aim of deepening our understanding of the process of 
cluster formation and evolution. 

A3921 is classified as a regular cluster of richness 2 
(Abell, Corwin & Olowin ITMty at redshift z = 0.094. It 
was detected first in X-ray by the Al survey of HEAO-1 
(Kolvalski et al. ! !)S Q) and was subsequently observed in 
pointed observations by the Einstein, ROSAT and GINGA 
satellites. The latter two observations have been analysed 
by Arnaud et al. (1996) who showed that A3921 con- 
sists of two merging components with a mass ratio of 1 
to 3. Their analysis indicated that the smaller subclus- 
ter is located at a distance of 0.8 hj^ 1 Mpc to the West 
of the main cluster. The ROSAT/GINGA observation in- 
dicated that the region between the subcluster and the 
main cluster has a higher temperature, however, the poor 
spatial resolution and limited energy coverage of ROSAT 
prevented deeper interpretation. A3921 was also observed 
with ASCA. White (2000) measured a global tempera- 
ture of 5.73lg'23 keV and a strong cooling flow of mass 
deposition rate M = 319 M Q yr" 1 . The temperature pro- 
file obtained by White (2000) is fairly isothermal, with 
an indication of a higher temperature at the edge of the 
cluster. However, the spatial resolution of the satellite did 
not allow a precise analysis. All published indications thus 
point to a merging cluster in an intermediate phase of the 
merger, just before the core passage. It was on the basis of 
this information that A3921 was selected for our sample 
of merging clusters of galaxies. 

In this paper we use the high sensitivity and good spa- 
tial resolution of XMM-Newton to gain new insights into 
this system. Our conclusions are the result of a comparison 
of our X-ray observation with numerical simulations and 
the new multi-object spectroscopy and VRI-band imag- 
ing data obtained at the ESO telescopes and detailed in 



an accompanying paper (Ferrari et al. 2004; hereafter 
F04) , and allow us to build an entirely new picture of the 
dynamical state of the system. 

Throughout the paper we assume Hq = 75 km s -1 
Mpc -1 , f2 m = 0.3, f^A = 0.7. In this cosmology, at the 
redshift of the cluster, 1 arcmin corresponds to 98 kpc . 

2. Observations and data preparation 

2.1. Observations 

A3921 was observed for 30 ks in October 2000 by XMM- 
Newton as a Guaranteed Time target. In this paper only 
data from the European Photon Imaging Camera (EPIC; 
Struder et al. 120011 Turner et al. 2001) are considered. 
Calibrated event list files were provided by the XMM- 
Newton SOC the 19 of June 2001. All observations were 
obtained with the MEDIUM filter in Full Frame mode. 
Throughout this analysis single pixel events for the pn 
data (PATTERN 0) were selected, while for the MOS data 
sets the PATTERNS 0-12 were used. 

In order to estimate the background level and its vari- 
ability with time because of flares, we produced light 
curves in 100s bins in a high energy band (10.0-12.0 keV 
for MOS and 12.0-14.-0 keV for pn). The mean level of 
the light curve was very regular, thus there are no flares 
and for the scientific analysis, we used the total exposure 
time of the observation, which was 30073.4 s, 30219.1 s, 
and 22398.1 s for MOS1, MOS2, and PN, respectively. 

The background estimate was obtained using the 
blank-sky event lists described in Lumb et al. l|2002i . The 
blank-sky background events were selected using the same 
selection criteria (such as PATTERN, FLAG, etc.) used 
for the observation events and transformed in order to 
have the same sky coordinates (EPOCH J2000) as for 
A3921. This procedure ensures that the source and back- 
ground products come from the same region of the de- 
tector, reducing errors induced by the detector position 
dependence. By generating the light curve of the blank- 
sky background in the high energy (10-12 keV) band, we 
rejected time intervals displaying more than 15 counts per 
100 s bin in order to have the same statistical flare- induced 
background as A3921. 

The source and background events were corrected for 
vignetting using the SAS task evigweight, which imple- 
ments the method outlined in Briel & Henrv IT994|I . This 
allows us to use the on-axis response matrices and effective 
areas. 

To take into account the variation of the particle in- 
duced background, which dominates at high energy (> 5 
keV) and induces fluorescence lines (Al, Si, Cu, Au) from 
the shielding of the camera and the detector itself, we 
computed normalisation factors by considering the pho- 
ton (particle) statistics in the 10-12 keV (MOS) and 12- 
14 keV band (pn) between the source and the blank-sky 
background observations. The normalisations used in this 
work are 0.99, 1.00, 1.12 for MOS1, MOS2 and pn respec- 
tively. 
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Fig. 1. (a) EPIC counts image in the energy band 0.3-10.0 keV; the image is not corrected for exposure, (b) Iso- 
intensity contours of the EMOS1+EMOS2 adaptively smoothed image in the energy band 0.3-10.0 keV overlaid on 
the deep I band optical image obtained at the MPG/ESO 2.2 m telescope (see IF04|l in la Silla. The brightest galaxies 
are marked. Contours come from a non-background subtracted, non-exposure-corrected image, and are logarithmically 
spaced. The lowest contour level is at 0.6 (MOS1+MOS2) counts/pixel. Epoch is J2000. 



3. The cluster X-ray morphology 

3.1. X-ray image 



Figure 1(a) shows the EPIC image in the 30' field of 
view (FoV) of XMM-Newton obtained by combining the 
MOS1, MOS2 and pn images in the energy range 0.3- 
10.0 keV. The image is not corrected for exposure. The 
cluster displays a slightly elliptical shape, with an exten- 
sion to the west-north- west (WNW). The extension dis- 
plays two extended emission peaks which correspond to 
two of the three brightest galaxies in optical, as shown 
in Fig. 1(b) Here the X-ray iso- intensity contours of the 
XMM-Newton MOS adaptively smoothed image are su- 
perimposed on the deep I-band image obtained with the 
WFI instrument at the MPG/ESO 2.2 m telescope |F04l 
for details). The two galaxies are at 7' (690h^ kpc; RA 
= 22 /l 48 m 49 s , Dec = -64°23'10" (J2000)) to the north- 
west (BG2) and at 8' (800/ifg 1 kpc; RA = 22 ,l 49'"04 s , 
Dec = -64°20'35" (J2000)) to the WNW (BG3) from 
the central brightest galaxy BG1 (RA = 22 /l 49 m 58 s ; 
Dec=-64°25'46" (J2000)). The peak of the X-ray emis- 
sion is shifted by 17" to the NW with respect to the opti- 
cal position of BG1. This is larger than any attitude error 
of XMM-Newton, and note that the X-ray and optical 
emission of BG3 coincide well within 17", although the 
X-ray emission around BG3 is extended. The peak of the 
BG2 X-ray emission is also shifted, this time towards the 



SE, by 26". The diffuse X-ray emission detected around 
galaxy BG3 corresponds to the emission of the subcluster 
detected by Arnaud et al. I|1996|> with the ROSAT data. 
We also notice an interesting depression in the surface 
brightness, between BG2 and BG3. 

Several point and extended sources are also detected 
in the field of A3921. Some of them are the X-ray coun- 
terparts of galaxies in the cluster, others are unrelated 
objects. All point sources with flux greater than 10 -14 
ergs s" 1 arcmin -2 were masked throughout this analysis 
by checking the source detection file in the pipeline prod- 
ucts. Such point sources were masked with a typical radius 
of 30". Some sources marginally detected by the SAS de- 
tection algorithm were also masked. 

3.2. 2D j3 model fitting 

For this analysis we have followed the method applied 
to the merging cluster A1750, described in Belsole et 
al. (2003J). We fit the surface brightness distribution 
of A3921 at low energy (where the distribution is less 
temperature-dependent) with a (3 model, and quantify the 
deviation from this model. We use only the MOS camera 
because the large gaps in the pn camera are difficult to 
take into account, and could produce spurious structure. 

In order for the model to be a good representation 
of the main cluster, we masked the region connected 
to the north-western extension with a rectangle of size 
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Table 1. Best- fitting results of the 2D (3 model. 



Parameter 


best-fitting 


r c i (arcmin) 


2.80 


r C 2 (arcmin) 


2.13 


P 


0.77 


RA (J2000) 


22 h 50 m 00 s 


Dec (J2000) 


65°25'43" 


PA (rad) 


0.19 



10.3 x 13.0 arcmin, together with the mask used to excise 
point sources. The centre of the cluster was also masked 
with a circle of 1.'36 (the results are identical whether or 
not the central region is excluded from the fit). The best- 
fitting 2D-model matches the elliptical shape and orienta- 
tion angle of the main cluster. The best-fitting parameters 
are listed in Tabled they are in very good agreement with 
results found by Arnaud et al. |1996). 




Fig. 2. EPIC/MOS Gaussian-filtered image with super- 
imposed contours of the residuals after the subtraction of 
a 2D /3-model- see text for details. The residuals are traced 
at 3,5,7,10,15 and 20 er significance and the image is in a 
logarithmic scale. Epoch J2000.[see the electronic version 
of the paper for a colour figure] 

We quantified the significance of possible excess flux 
by subtracting the model from the image. The residuals 
are shown in Fig. [3] We detect extended residuals, roughly 
centred on the galaxy BG3, at more than 5tr significance, 
supporting the previous indications of the double nature 
of A3921. In this framework, BG3 would be the central 
galaxy of a subcluster which is already highly disrupted 
by interaction with the main cluster, in agreement with 
the analysis of Arnaud et al. I|1996ll . Within this extended 
residual structure there are two peaks, other than BG3, 
one of which coincides with the position of galaxy BG2. 



The residuals are elongated in the direction of the main 
cluster to the East and display a sharp edge towards the 
NW. We also detect excess emission in the centre of the 
main cluster, reminiscent of a cooling flow. This residual 
emission displays a fairly circular shape with an elongation 
to the east and compression of the isophotes to the west, 
in the direction of the large extended residual structure. 

3.3. ID analysis 

We also performed a one-dimensional analysis by extract- 
ing a radial profile of the main cluster, binning events in 
circular annuli centred on the emission peak. We excluded 
the side of the cluster corresponding to the residuals shown 
in Figure (we use the main region described below and 
illustrated in Fig. 03). In view of the central excess in the 
2D analysis above, it is obvious that a single /3-model will 
not be a good description of this radial profile. It was thus 
fitted with a double isothermal /3-model (the BB model, 
see Pratt & Arnaud 2002 ) . The model was convolved with 
the instrument response (Ghizzardi 2001 ) and binned into 
the same bins as the observed profile. The best-fitting 
values are: r c = 2.43±£^ arcmin, [3 = 0.77 ± 0.02, 
''cut = 1-581q'j3 arcmm j an d r c , in = 022 ± 0.07 arcmin 
(90 per cent errors for one interesting parameter), which 
gives x 2 = 223.2 for 179 degrees of freedom. It can be 
seen that the external f3 and r c values are in excellent 
agreement with the 2D analysis. 



4. Thermal structure 

In this Section, we will deal with the temperature struc- 
ture of the cluster. We will first investigate the global 
spectral properties of the cluster. Then, after calculating 
a temperature map, we will extract spectra in distinct re- 
gions as an independent check of the significance of the 
features we find. 

In Appendix A, we discuss the role of possible sys- 
tematic effects such as variations in abundance or column 
density, and differences in plasma codes. Our conclusion 
from these tests is that the temperature structure we find 
is robust to these sources of systematic error. 

The A^h value used for the temperature map compu- 
tation, and for any spectrum described in the following, 
is obtained by simultaneously fitting the MOS spectra 
only (between 0.3 and 10.0 keV) with an absorbed mekal 
model. This gives N H = 2.19i^4 2 x lo2 ° cm ~ 2 ( 90 P er 
cent errors for one significant parameter), which is ~ 25 
per cent lower than the Galactic value given by Dickey 
& Lockman (1990), but in excellent agreement with the 
PSPC result of N K = (2.2 ± 0.2) x 10 20 cm" 2 (Arnaud 
et al. H99ri|l . Possible systematics related to Nh variations 
are described in Appendix A. 
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4.1. Global spectrum 

In this Section we are interested in evaluating the global 
properties of the cluster, for comparison with previous 
works, and to quantify the effect of the merger on them. 

The spectral analysis discussed in this Sect, and the 
following is performed by extracting spectra from the 
weighted source and background (blank-sky) event lists. 
Spectra were extracted in each region, and the background 
spectrum was normalised and subtracted from the source 
spectrum. The spectra of the three cameras were fitted si- 
multaneously in the (0.3/0.5 - 10.0) keV range (MOS/pn). 

We extracted a spectrum in a large circle of radius 
11', corresponding to ~ 0.65 r2oo • We will refer to this 
spectrum as the global spectrum. 

After considerations about possible systematics errors 
due to different plasma codes and/or specific element 
variations (tests are described in detail in Appendix A), 
we fitted the global spectrum with a single and a two- 
temperature mekal model. Results are listed in Table [5] 
The single-temperature model gives a reduced x 2 — 1.36, 
clearly not satisfactory. Adopting this model, the global 
temperature, kT = 4.34±0.10 keV (90 pe cent errors), is in 
agreement with the GINGA (kT = 4.9 ±0.6) and ROSAT 
(kT = 4.1 ±1.0) results of Arnaud et al. ©96), but is sig- 
nificantly lower than the ASCA results (kT = 5.73 ± 0.24; 
White 2000). 

With the single temperature fit to the global spectrum, 
we observe large residuals corresponding to the Fe L mul- 
tiplct (see the top panel of Fig. |3J), reflecting the exis- 
tence of multi-temperature gas as observed in the tem- 
perature map above. The largest residuals are centred at 
~ 0.8 keV. If the Fe-L-shell lines are used as a thermome- 
ter (Belsole et al. 2001, Bohringer et al. 2002), we might 
expect to observe at least another component at ~ 0.6 
keV. We thus fitted the spectrum with two-temperature 
mekal model, obtaining significantly better fit (reduced 
X 2 =l.l). We observe that both the dominant tempera- 
ture and the chemical abundances are higher in this case 
(but see Appendix A). The two-temperature fit is shown 
in the bottom panel of Fig. 13 The low temperature (kT = 
0.65 keV) and luminosity (L x = (1.7 ±0.2) x 10 43 erg s" 1 , 
more than an order of magnitude lower than the first com- 
ponent) of the low-temperature thermal component may 
originate from the central galaxy. To test this hypothe- 
sis, we fitted the central bin of the temperature profile 
(see Sect. 14.4]) . e.g. the region within 30 arcsec, with one 
and two mekal models. The spectrum is well fitted with 
a 1-temperature model of kT = 5.44 keV (not requiring 
a second component), indicating that any galaxy compo- 
nent is negligible with respect to the cluster emission. The 
low-temperature component more likely arises from the 
fact that we are fitting a multi-temperature plasma with 
a two-temperature model (see also Appendix). 




1 Here the virial radius is defined as r2oo = h(z) 1 x 3.69 x 
(T/10 keV) 172 ^ 1 Mpc, where h{z) = (O m x (1 + zf + 
A) 172 and the relation is scaled from Evrard, Metzner & 
Navarro fHMt . 



channel energy (keV) 

Fig. 3. Global spectrum integrated in a circular region 
of radius 11.'2, corresponding to ~ 0.65 r 2 oo- The pn 
is in grey the MOS in black. Top: The folded model is 
the single-temperature mekal. One can observe the large 
residuals at ~ 0.8 keV. Bottom: The folded model is the 
two-temperature mekal. The spectra are binned to have 
a S/N ratio of 5<r above the background. 



White (2000) detected a cooling flow with a large mass 
deposition rate, and we see a residual in the core in the 2D 
/3-modcl fit. However, our temperature map (and profile) 
clearly shows that there is no cool gas in the core. 

4.2. Temperature map 

The temperature map was obtained by applying the multi- 
scale spectro-imaging algorithm described in Bourdin et 
al. (2004). Details of the application of this technique to 
XMM-Newton data can be found in Belsole et al. ( 2004 ) . 
Although we did not include the pn camera because of 
technical considerations, in the case of A3921 we have 
more counts than we had for A1750 and the temperature 
structure is mapped in finer detail. To enhance the stabil- 
ity of our results we impose a (conservative) lower limit 
of 1000 photons per pixel at each spatial scale. Both the 
particle and X-ray backgrounds are modelled in this ap- 
proach, allowing us to achieve typical errors of ~ 0.1 keV 
per pixel (see Bourdin et al. 120041 for details). The fit in 
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Fig. 4. Temperature map obtained by applying the wavelet algorithm described in Bourdin et al. |2004). The contours 
are obtained from the MOS adaptively smoothed image in the energy band 0.3-10 keV, and are logarithmically spaced, 
as in Fig.H The first contour corresponds to 0.6 (MOS1+MOS2) counts/pixel. Image coordinates epoch is J2000. [see 
the electronic version of the paper for a colour figure]. 



each pixel and spatial scale was performed with the metal- 
licity fixed to the best-fitting single-temperature MEKAL 
fit to the global spectrum, and the column density fixed 
to the best-fitting value derived from the MOS spectral 
fit to the main spectrum (Ah — 2.19 x 10 20 cm -2 ; see 
Appendix A for details). 

The temperature map obtained from this wavelet anal- 
ysis is shown in Fig. 01 The overlaid contours are the same 
as in Fig. 1(b) The overall appearance of the temperature 
structure is strongly asymmetric. We observe an extended 
hot region which appears related to the extended residuals 
seen in Fig. The core of the main cluster is also hotter 
than the surroundings, at least to its south-western side. 
The hottest region of the map is oriented nearly parallel to 
the line joining the centre of the main cluster and a point 



between the two dominant galaxies to the west, and in- 
deed terminates in a region between BG2 and BG3. In the 
direction of the main cluster, at 1.5' (~ 150h^ kpc) from 
the centre towards the west, the hot emission is more el- 
liptically shaped and follows the X-ray emission isophotes. 
The eastern side of the main cluster is fairly isothermal, 
with an average temperature of 4.9 keV. The overall tem- 
perature structure is difficult to understand if the main 
cluster and the secondary structure are just starting to 
collide. 

We also detect a hot region to the south of the core, 
at around 5.3' (RA = 22 /l 49 m 35.8 s , Dec = -64°30'25'.'35 
(J2000)). We did not find any correspondence between this 
hot spot and any X-ray emission peak. We notice that the 
emissivity contours are distorted at this location. In an 
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Fig. 5. Figure illustrating discrete regions from which 
spectra were extracted. The regions were defined on the 
basis of the residual morphological and temperature struc- 
ture [see the electronic version of the paper for a colour 
figure]. 

attempt to quantify these emissivity distortions, we gen- 
erated a radial profile where we azimuthally integrated 
counts within a sector (of angle 60°) centred at the clus- 
ter centre and including the hot region. We did not find 
significant deviation from a /3-model. We also compared 
with the deep optical image l|Ff)4J> . but there is no opti- 
cal counterpart which could explain a hot emission at this 
location. Thus everything seems to indicate that this hot 
emission is related to the merger event. 

4.3. Region analysis 

We now turn to an assessment of the significance of the 
temperature structure found in the temperature map. In 
the following analysis, we will concentrate on spectral fits 
to discrete regions with mekal models. Except where 
noted, the column density is fixed to the best-fitting value 
derived from the MOS spectral fit to the main spectrum 
(Ah = 2.19 x 10 20 cm~ 2 ). The individual regions are illus- 
trated in Fig. 

4.3.1. Discrete region 1: the main cluster 

We extracted a spectrum for the main cluster only by 
selecting a sector between PA = 71° and 268° (East to 
North, illustrated in Fig. EJ, e.g. excluding the western 
extension and hot emission bar. We will refer to this spec- 
trum as the main spectrum. We fitted the main spectrum 
with single and double MEKAL models, for which the re- 
sults are detailed in Table [3 The main region has a tem- 
perature which is slightly lower than the global tempera- 



ture, a direct consequence of excluding the hot bar (see 
Fig. 0}. This temperature is also lower than the average 
temperature on the eastern side of the temperature map. 
Part of this discrepancy can be explained by the fact that 
we did not use the pn camera when deriving the tempera- 
ture map. Once again the two-temperature model is a bet- 
ter representation of the data (F-test probability > 99.9 
per cent), suggesting the existence of multi-phase gas even 
on the side of the cluster which looks the most dynamically 
relaxed. The ratio between the bolometric X-ray luminos- 
ity found in the global and main regions is L 9L \ ohal / U^ am 
= 2.0 which is comparable to the ratio of the two areas of 
integration. 

4.3.2. Discrete region 2: the hot bar 

We next extracted the spectrum of the hot bar observed 
in the temperature map. The region is an ellipse of semi- 
minor axis 98", semi- major axis 174", and position angle 
45°, centred on RA = 22 /l 49 m 14 s 813, Dec = -64° 23. '49'.' 51 
(see Fig. El . Results are listed in TableEl As expected, the 
temperature of the hot bar is significantly higher than the 
global and main cluster regions. Possible systematics can 
arise because of absorption variations (see Appendix A). 
To test this effect, we left the iVn free to vary, obtaining a 
best-fitting Nn lower than (but consistent with) the aver- 
age value and an even higher temperature (although with 
larger errors). Thus if there are any systematics related 
to absorption variations, the approach of fixing the Nn 
for the derivation of the temperature map is somewhat 
conservative. 

When a second thermal component is added to the 
model, the dominant component reaches a temperature of 
7.25 keV. The F-test indicates that the two-temperature 
model is a better representation of the data, with a proba- 
bility of 3.14 x 1CF 3 that the second component is detected 
by chance. FigureQJhas already graphically illustrated the 
multi-temperature nature of the hot bar. 



4.3.3. Discrete region 3: the western region 

In order to further verify the asymmetric temperature dis- 
tribution represented by the hot bar, we extracted a spec- 
trum in the region to the west, beyond the hot bar, as 
shown in Fig. We excluded the X-ray emission coinci- 
dent with the two brightest galaxies BG2 and BG3. 

This region is slightly cooler than the main region (but 
consistent within the errors) . Once again, we observe that 
a two-temperature thermal model is a better representa- 
tion of the data than a single temperature model (F-test 
probability >99.9 per cent). We also observe an increase 
in relative abundances, which might suggest chemical con- 
tribution from a more metal-rich object at the location of 
the residual structure. 
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Table 2. Best-fitting results of discrete regions. 



boi 
X 



( 2 /d.O,f. 



region 



Tl 

(keV) 



Zl 
solar 



Normf 



j+U.UV 
-0.07 



26 +u - ua 

u - ZD -0.03 
n QQ+0-04 
U-O»_o .04 

u - ZJ -0.04 

o.39±g;gS 
0.28lg;ii 

0.36l°£ 
0-32±°;i| 
0.731°;^ 



1 - oo -0.02 

i 7n+ - 04 

Q7+ 01 
u - y ' -0.01 
7+0.02 
-0.02 
-,+0.04 
-0.04 
,+0.03 
-0.04 
j+0.47 
5 -0.46 
-+0.63 
-0.34 



T2 

(keV) 



Z2 
solar 



Norm2 

n 



(10 44 erg s" 1 ) 



°- Z4 -o.io 

K =4+0.10 

2.65l° ;l° 
3.80l° ;i° 
n 4n +002 

u - w -0.02 

n 42+ - 05 

u -* z -0.04 

n 2n+ ' 02 

u - zu -0.02 
n 99+O.O3 



global 
global 
main 
main 
hot bar 
hot bar 
western 
western 



4.38^ 

c 18 +0.18 
O.18_ 17 

4.i6±S:lS 

4.90t° : * 
6.30l°;^ 

7 9 r+0.79 
' - zd -0.62 
o fin+0 39 
a.DU_ 29 



4.59 



+0.61 
0.48 



0. 871 
1.20j 

1. u] 

7.( 

6.35^ 



0.65 



+0.04 
0.04 



0.17 



+0.15 
0.06 



j+0.05 
'-0.06 



0.14 



+0.17 
0.06 



0.10] 



x 10" 



0.64 



+0.22 

0.17 



0.67 



+0.09 
0.01 



I +2.96 



3.25t};2? x 10" 



1551.1/1140 
1256.1/1136 
1190.2/898 
1063.1/894 
708.2/707 
687.0/704 
419.9/384 
387.6/381 



Errors are quoted at 90 per cent confidence for one interesting parameter. For the hot bar and the western regions the 
abundances of the two thermal components were tied.* The normalisation is in units of 10- 14 /(47r(D A x (1 + z)) 2 ) cm , 
where Da is the angular size distance to the source (cm), as defined in the mekal model, in XSPEC. 



4.4. Temperature profile 

The gas density and temperature profiles are key observa- 
tional quantities for the derivation of many other cluster 
characteristics. In particular, these quantities are essen- 
tial for the derivation of mass and entropy profiles. In 
this Section, we will calculate the temperature profile of 
A3921, and assess possible systematics. Most analyses of 
more relaxed clusters assume single temperature models, 
and so we ignore multi-temperature effects in the present 
analysis. 

We first obtained a temperature profile in the main 
region, (Fig. |SJ. The centre of the profile was fixed to the 
peak of emission of the main cluster. The radial tempera- 
ture bins were defined by requiring a signal-to-noise of 5<r 
in the 2.5-6.5 keV radial profile, with a lower limit of 30" 
to minimise Point Spread Function (PSF) effects. (The 
annuli thus have widths greater than or equal to the di- 
ameter enclosing 70 per cent of the energy for the on-axis 
PSF). We were able to extract spectra in 10 annuli, out to 
a distance which corresponds to ~ 0.5 r2oo- We fitted each 
spectrum with a single mekal model where the column 
density was fixed to the value used for the temperature 
map and the discrete regions (JVh = 2.19 x 10 20 cm" 2 ). 
Chemical abundances were left as free parameters out to 
annulus 8 and then fixed to the fitted value obtained for 
annulus 5. 2 The three camera spectra were fitted simulta- 
neously, ignoring channels below 0.3 for MOS and below 
0.5 for pn. We found excellent \ 2 f° r a U our bins, indi- 
cating that a second thermal component is not required if 
the spectrum extraction region is sufficiently small. The 
resulting projected temperature profile is shown in Fig. EI 

The quiescent background level of XMM-Newton can 
vary by ±10 per cent. We thus calculated temperature 
profiles with the background at plus and minus 10 per 
cent of our adopted normalisation. The effect is negligi- 
ble, both of these profiles are consistent with the nominal 
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When the abundances were fixed to 0.25, the fitted tem- 
perature were statistically consistent with the values obtained 
when the abundances were left free. 



2 4 6 8 

Radius ( arcmin ) 

Fig. 6. Projected radial temperature profile of the main 
and global regions. The lines ovcrplotted on the main re- 
gion profiles are the best-fitting polytropic models. The 
values of To and 7 for all of these models agree within 
their la errors (see text for details). 



profile at la, as can be seen in Fig. El For easier compar- 
ison with previous work, we also extracted a temperature 
profile from annular regions encompassing the entire clus- 
ter (the global region). This is also shown in Fig. [HJ 

There is no sign of a cool core region, and all the pro- 
files drop almost monotonically to about 50 per cent of 
the central value at the detection limit. The global region 
profile is slightly hotter and flatter than that of the main 
region, as expected since this region includes the emis- 
sion from the hot bar. However, the temperature still de- 
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Fig. 7. Integrated total mass profile (points with errors), 
using data from the unperturbed part of the cluster. The 
solid line is the best-fitting NFW profile (xt ~ 4.5). 



clines with radius. This is different from the flat tempera- 
ture profile found by White (2000) with ASCA. We have 
fitted a polytropic model (e.g., Markevitch et al. 0,998), 
T = Tq x p 7_1 , to the main region profiles. The gas den- 
sity profile was described by the BB model discussed in 
Sect. 13.31 with all parameters fixed. We find that the nom- 
inal main region profile has To = 6.07 ±0.03 keV and poly- 
tropic index 7 = 1.16 ± 0.01. The values of To and 7 for 
the main region profiles with background at plus/minus 
10 percent of nominal agree within their la errors. 

4.5. Luminosity 

Numerical simulations (Ritchie & Thomas. 2002 Randall, 
Sarazin & Ricker I2002fl suggest that the effect of the 
merger is to boost the luminosity and the temperature 
for a short time, and this can have an obvious effect on 
scaling relations which involve these quantities. Arnaud & 
Evrard Ijl999|) found Lx oc T 2 88 using a sample of clus- 
ters with weak or absent cooling flow signatures. Their 
relation has an intrinsic dispersion of 0.13 in log Lx at a 
given temperature. Taking the simple onc-tcmpcraturc fit 
to the global and main regions in our analysis and scal- 
ing the Lx - T relation of Arnaud & Evrard 1|1999[) for 
our cosmology, we find that A3921 is a factor of two too 
luminous for its temperature. We further note that this 
cluster is an outlier in their analysis (see their Fig. 1), but 
less significantly so because of the larger errors on their 
temperature measurement. 



5. Mass analysis 

Combining a gas density (Sect. I3.3fl and a temperature 
(Sect. |OI) profile, we can derive a total gravitational mass 
profile under the assumptions of hydrostatic equilibrium 
and spherical symmetry. Given the observed ellipticity 
and substructure (Sect. I3.2|l and the strong merger sig- 
natures (Sect. l4~2*)l . it is very likely that the cluster is not 
in hydrostatic equilibrium. We have thus tried to min- 
imise these effects by excluding the regions most strongly 
affected by the merger event. This approach has often 
been used for clusters with substructure (e.g., Henry, Bricl 
& Nulsen I1993|l . and the approach has been shown to 
give reasonable results when applied to simulated clusters 
rSchindler lTTMl) . 

We thus use the main region to derive the mass profile. 
The mass was calculated at each radius of the temperature 
profile with the Monte Carlo method described in Pratt & 
Arnaud ( 2003 ) , which takes as input the parametric model 
for the gas density and the measured temperature profile 
with errors. A random temperature is calculated at each 
radius of the measured temperature profile, assuming a 
Gaussian distribution with sigma equal to the la error, 
and a cublic spline interpolation is used to compute the 
derivative. Only profiles corresponding to a monotonically 
increasing mass gradient are kept: 1000 such profiles were 
calculated. The final gravitational mass profile is shown 
in Fig. □ 

For this cluster at least, differences in temperature pro- 
files have a negligible effect on the resulting mass profile. 
Mass profiles calculated using all of the temperature pro- 
files detailed above in Sect . 14.41 (Including that of the global 
region which contains all the disturbed temperature struc- 
ture) all agree within the la errors. 

We have tried to fit the mass profile with several dif- 
ferent functional forms for the dark matter density. The 
Moore et al. ( 1999) and King approximation to an isother- 
mal sphere models are unconstrained. An NFW model 
(shown in Fig.0) yields c = 4.3 and r s = 337hj^ kpc, but 
has a reduced xl = 4.5. Thus either i) the density profile 
of this cluster is not well described by the NFW model, or 
ii) even in the relaxed-looking main region, the gas is not 
in hydrostatic equilibrium. In any case, the poor fit argues 
in favour of the cluster being dynamically perturbed. 

It is important to understand the effect of mergers on 
mass calculations. The NFW fit is particularly poor in 
the inner regions (r < 150kpc), but beyond this radius 
it is a rather good description of the X-ray mass profile. 
This suggests that it might be interesting to compare the 
total optical and X-ray masses within the fiducial opti- 
cal radius of 0.34/ifg 1 Mpc i|F04jl . The NFW fit yields 
9.8 x lO 13 /^ M Q interior to this radius, while a linear 
extrapolation of the X-ray mass profile in the log r-log 
M plane gives (9.9 ± 1.2) x 10 13 h^M G . The optically- 
derived mass of the main cluster within this radius is 
2.1±q% x 10 u h^M e . Thus at this radius the different 
mass estimators agree within their respective errors. 
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6. The centre of the main cluster 

In Sect. 13.21 we found that the centre of the cluster shows 
an excess of emission over a /3-model, and that this re- 
gion also has a relatively high temperature (5.4 keV), as 
shown by the temperature profile and map, which also 
suggest that this central region is in a multi-phase state. 
In principle, while the high temperature of the first bin in 
the temperature profile might be produced by an active 
galactic nucleus (AGN) in the central galaxy BG1, we can 
find no evidence for a second component in the spectrum 
from this bin. In addition, the large extent of the spa- 
tial residuals (2 arcmin ~ 130 kpc) make it unlikely that 
what we observe could be produced by BG1 alone. We 
thus conclude that the high central temperature is likely 
a consequence of the merging event. 

We looked at the central 450 kpc (a 4f5 x 4.'5 square 
region) in more detail by extracting images in 5 different 
energy bands. This allows us to compare simultaneously 
the spectral and spatial structures. In Fig.|B|we show from 
the top left to the bottom centre the images in the (0.2- 
0.7) keV, (0.75-1.25) keV, (1.3-1.8) keV, (2.0-4.5) keV and 
(4.7-8.0) keV bands obtained with the MOS1 camera. The 
images have been smoothed with a Gaussian filter of a = 
11". The cross is centred on the emission peak in the (0.2- 
0.7) keV image. In order to check for PSF effects, we also 
produced images in the same energy bands for a bright 
point source approximately 9' off-axis in the same field 
(see Fig. bottom right) - note that the X-ray peak of 
A3921 is 4' off-axis. The point source image is shown in 
the insert for each energy band. 

Despite the fact that the point source is further off-axis 
than the X-ray peak of A3921, we note that its centroid 
position does not change with energy (except perhaps for 
the hardest energy band, where we measure a small shift 
of 3'.'3 toward the East). The centroid of the A3921 emis- 
sion peak changes significantly with energy. The largest 
shift (excluding the hardest energy band) is 15" for the 
(1.3-1.8) keV image. We also observe isophotal orientation 
variation. 

The lowest-energy image is the most sensitive to small- 
scale changes in the column density, and the lower than 
Galactic Nh we have found in the global spectral fit sup- 
ports this hypothesis. This difference in Nh values is a 
possible explanation for the shifted centre position with 
respect to other energy bands. The (0.75-1.25) keV im- 
age covers to the Fe L multiplet. The centroid shift and 
the North-South elongation of the emission in this energy 
band is possibly connected to abundance variations, or, 
more likely, temperature variations. The shift in the hard- 
est energy band is partially due to the PSF, as observed 
for the point source (in the insert). However, PSF effects 
cannot explain the whole shift. For the other two energy 
intervals, we cannot find an obvious physical explanation 
for the centroid shift, apart from perturbations due to the 
merger event. The Chandra observation of A3921, to be 
performed shortly, will shed more light on this issue. 




Fig. 8. Images of the central 200 kpc of A3921. All images 
are smoothed with a Gaussian filter of a = 11". From top 
left to bottom centre the energy bands are: 0.2 - 0.7 keV; 
0.75 - 1.25 keV; 1.3 -1.8 keV; 2.0 -4.5 keV and 4.7-8.0 keV. 
The insert in each figure shows the corresponding point 
source image. The bottom right image shows the regions 
used for the image extraction: in black for the point source 
and in white for A3921.[see the electronic version of the 
paper for a colour figure] 



7. The bright galaxies BG2 and BG3 in X-ray 

We have also examined the X-ray properties of the bright 
galaxies BG2 and BG3. We extracted the spectrum of each 
galaxy in circular regions of radius of 20" and 16" for 
BG2 and BG3 respectively. As a background we selected 
an annular region of inner radius 33" and outer radius 
53" which takes into account the cluster emission at the 
location of the galaxies. 

There are 300 net MOS+pn counts from BG2 between 
0.3 and 5.0 keV, severely limiting the constraints we can 
put on the spectral parameters. We fitted simultaneously 
the three camera spectra with an absorbed mekal model. 
We fixed the galactic absorption to the global X-ray value 
(2.19xl0 20 cm -2 ), and let the temperature and abun- 
dance vary. The best-fitting gives kT = l.lOig'is anc ^ a 
low mctallicity of 0.05(< 0.14) solar, with a \ 2 = 12.5/21 
d.o.f (errors are quoted at la confidence). If we fix the 
metallicity to be 0.3, the % 2 rises to 15.6 for 22 d.o.f, and 
the best-fitting temperature is higher (1.34io - 23 keV). 

We are also interested if the X-ray emission has a 
non-thermal origin and thus we also fitted the spec- 
trum with a power-law model. We obtain a photon in- 
dex 7 = 2.26l5'26' consistent with X-ray results of ac- 
tive nuclei (e.g. Brinkmann, Papadakis & Ferrero E004). 
The x 2 is 11.1/19 d.o.f. In this case, the thermal compo- 
nent best-fitting parameters are kT = 0.9(±0.4) keV and 
Z — 0.7(> 0.13). The thermal and non-thermal compo- 
nent have similar normalisations. 

The thermal and non-thermal models fit the data 
equally well, but it has to be admitted that the data are 
not very constraining. However, the non-thermal model 
implies that BG2 harbours an active nucleus. A3921 was 
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observed with the Australia Telescope Compact Array 
(ATCA) at 2.5 GHz and 1.5 GHz. The 2.5 GHz radio map 
(M. Hardcastle, private communication) shows faint radio 
emission in the centre of the optical counterpart of BG2. 
No radio emission is detected in the lower frequency obser- 
vation (1.5 GHz), suggesting that the source may indeed 
be an AGN. We conclude that the X-ray emission from 
BG2 may be non-thermal in origin, more likely coming 
from a central active nucleus. The weakness of the X-ray 
and radio emission could be a consequence of the merging 
event. If the gas surrounding the core of BG2 is stripped 
by the merger, the active nucleus would be starved of fuel. 

The X-ray spectrum of BG3, on the other hand, looks 
more thermal. The statistical quality here is slightly bet- 
ter, with 406 MOS+pn photons in the (0.3-3.0) keV energy 
range. We fitted simultaneously the three camera spectra 
with an absorbed mekal model, the Nn fixed as above. 
The best-fitting temperature is kT = 1.07±g;|J keV (la- 
confidence), and the abundances Z = O.lllg'J? solar, for 
a x 2 = 25.5/26 d.o.f. With this model we calculate a 
bolometric luminosity of 6.0xl0 41 erg s _1 . Since contam- 
ination from the cluster may not have been completely 
accounted for by the background region we have cho- 
sen, we also fitted the spectrum with a two-temperature 
model where the high temperature was fixed to the cluster 
temperature at the position of BG3. We found that this 
model is a slightly (not significant) better representation 
of the data (x 2 = 24.8/25). The temperature and abun- 
dances are in agreement with the 1-temperature model, 
kT = 0.90±g;ff keV, Z = 0.48tgJt 

The relatively high X-ray luminosity suggests that 
BG3 belongs to the category of X-ray bright galaxies (e.g. 
Forman, Jones & Tucker 119851 Trinchieri, Fabbiano, & 
Canizares 1986), for which the diffuse X-ray emission from 
the IGM dominates over the emission of a population of 
non-resolved sources (mainly low mass X-Ray binaries; 
e.g., Sivakoff, Sarazin, Irwin 2003). However, the emission 
around BG3 looks more extended than that of a normal 
early-type galaxy. With our data we cannot exclude the 
possibility that the relatively high X-ray temperature is 
due to the emission from a small galaxy group with BG3 
at its centre, and which presumably has not yet collided 
with the main cluster or subcluster. It is difficult with 
the present data to separate any possible group emission 
from the extended emission of the X-ray residuals (see 
Fig. EJ). However F04 do not find other bright galaxies 
around BG3 in the projected redshift space. Deeper opti- 
cal observations in the region round BG3, together with 
the upcoming Chandra observation, will help shed light on 
the nature of BG3. 



8. Discussion 

8.1. A new picture 

The morphological analysis of the XMM-Newton obser- 
vation (Sect. l3~2*|) has shown unequivocally that there are 
two peaks in the X-ray emission from the system. The 



bulk of the X-ray emission is centred on the cD galaxy 
BG1. The residuals above a simple 2D j3- model fit to the 
main cluster are centred on BG3, which is clearly visible 



as an extended source even in the raw image (Fig. 1(a) I. 

The temperature analysis (Sect. l4~2*|l has shown a com- 
plicated temperature structure to the West of the main 
cluster, coincident with the residuals from the morpho- 
logical analysis. In particular, we detect a "bar" of signifi- 
cantly higher temperature gas (kT from 7 to 8 keV) which 
describes an angle which is approximately parallel to the 
line joining BG1 and BG2. If we draw a line in the SE-NW 
direction, passing through BG1, we see that the tempera- 
ture structure above the line is far more uniform than that 
below the line. Our analysis in discrete regions combined 
with the stable temperature fit after considering several 
possible systematic effects firmly demonstrates that the 
cluster temperature distribution is strongly asymmetric, 
with the temperature of the hot bar being more than 2 
keV higher than the less perturbed side of the cluster. 

The centre of the main cluster shows strong evidence 
for having been perturbed, both in shifts of the X-ray cen- 
troid with energy (Sect.EJ), and in the offset between the 
position of BG1 with that of the overall X-ray centroid. 
We note that the offset is in the same SE-NW direction as 
the hot bar. Furthermore, the monotonically-rising tem- 
perature profile towards the centre is neither in agreement 
with other XMM-Newton observations of relaxed clusters, 
nor with the concept of a cooling core region, a natural 
conclusion from the surface brightness analysis fSect. l3~2*l 
and l3.3|l alone. Moreover, the central regions of the mass 
profile (derived from the relaxed-looking main region) are 
not well fitted by the NFW profile, which is a reasonable 
representation of the mass distribution at larger radii. 

The optical observations of IF04I have given firm ev- 
idence for two peaks in the galaxy density distribution, 
one centred on BG1 (main component, A) and the other 
centred around BG2 (component B), the two components 
being at the same mean redshift. 

If we did not know the galaxy redshifts, it would be 
tempting to conclude that the residual emission we see is a 
secondary diffuse X-ray component with BG3 at the cen- 
tre, as was earlier suggested by Arnaud et al. 1)1996). BG3 
lies at a significantly different redshift from that of the two 
main galaxy peaks. From the current spectroscopic results 
of F04 there is no evidence at the moment for a concentra- 
tion of galaxies around BG3 in the redshift space. Other 
galaxies are measured to have the same cz as BG3, but 
they are located more in the main cluster central region 
(see IF04I for details). However, the X-ray temperature of 
BG3 is rather high (Sect. EJ), which could mean that it 
is associated with a group. Deeper optical observations 
should resolve the issue. 

8.2. Tentative interpretation 



Previous detailed work on this cluster (Arnaud et al. H996[) 
suggested that the system is in an early merger phase, with 
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the cluster B on its first infall towards the main cluster. 
With the present XMM-Newton observation and in con- 
junction with results presented in F04J we have several 
pieces of evidence that this interpretation is incorrect. 

The existence of the hot bar and particularly its orien- 
tation is the strongest evidence that we are not observing 
an early stage of the merger event. If we were, we would 
expect any compressed or shocked region to be orthogo- 
nal to the observed orientation, as suggested by numeri- 
cal simulations of merging clusters (Roettiger et al. 119971 
Ritchie & Thomas 0003 Ricker & Sarazin I2UUT1 Teyssier 
2002). In fact, we can compare our surface brightness dis- 
tribution and temperature map (Fig. ^ directly with the 
simulations by Ricker & Sarazin (|2001() and in particular 
with the bottom set of panels in their Fig. 7. These sim- 
ulations describe a 1:3 mass ratio off-axis merger with an 
impact parameter b of 5r s (where r s is the scale radius 
of the NFW profile describing the density distribution of 
the main cluster) . Our observed cluster surface brightness 
distribution resembles closely, but not exactly, the pre- 
dicted X-ray surface brightness distribution 1 Gyr after 
maximum luminosity. At the same time, our temperature 
map resembles their emission-weighted temperature map 
at maximum luminosity (0 Gyr). 

We can estimate the time since closest core passage 
by further comparison with the simulations of Ricker & 
Sarazin, taking the above mentioned 1:3 mass ratio, b = 
5r s merger case. In Sect 14.31 with the adoption of the 
Lx — T relation of Arnaud & Evrard l|1999[) for nearby 
clusters, we found that the cluster is roughly 2 times more 
luminous for its temperature. Let us consider the plot of 
X-ray luminosity versus time in the Ricker & Sarazin's 
Fig. 8, taking the quiescent luminosity to be the value 
4 Gyr after maximum. The simulated cluster is twice as 
luminous at a time ~ 0.4 Gyr after maximum luminosity. 
Note that the average emission-weighted temperature of 
the simulated cluster at this time has already dropped to 
approximately the quiescent temperature. 

An alternative estimate of the merger age can be ob- 
tained by assuming that the NW edge of the bar is the 
shock front associated with the off-axis collision between 
the main cluster and the sub-cluster associated with BG2. 
In the Ricker & Sarazin simulations, the shock forms when 
the outer regions of the the objects start to interact (their 
Fig. 6). Let us make the simple assumption that the shock 
front initially formed at the bottom (South) edge of the 
hot bar we observe today. Then we can speculate (assum- 
ing interaction in the plane of the sky) that the shock 
front has travelled a distance equivalent to the length of 
the hot bar, i.e., about 900 /175 _1 kpc. The pre-shock gas 
detected around BG2, in the northern outskirts of the hot 
bar, is at a temperature kTi=4.8 keV. The shocked gas 
at the NW edge of the hot bar is at kT 2 = 7.8 keV (here 
we adopt the value obtained via the temperature map), 
yielding a Mach number for the shock of M—1.6. From 
the Mach number and considering that the shock is mov- 
ing at the speed of the gas behind the shock itself, we 
calculate that the shock is moving at u ~ 1800 km s _1 . 



Combining this result with the length of the bar, we find 
that the shock front has travelled 0.9 Mpc at a speed 
of 1800 km a , implying a merger timescale of t < 0.5 
Gyr. This rough estimate is in very good agreement with 
the result discussed above. 

We thus suggest that we are observing an off-axis 
merger between two unequal mass objects, where subclus- 
ter A (the main cluster) represents the more massive com- 
ponent and subcluster B (around BG2), having come from 
somewhere in the SSE direction (where we observe the be- 
ginning of the hot bar), has already passed beyond closest 
core passage and is progressing on an outgoing trajectory 
towards the NW. This is a natural explanation for the ob- 
served surface brightness and temperature structure. Our 
interpretation is further bolstered by the lack of observ- 
able merger signatures in the radial velocity distribution, 
which would seem to suggest that the motion is taking 
place in the plane of the sky (see F04 for details). If sub- 
cluster B has already passed through the atmosphere of 
subcluster A, we would expect that the effect of the inter- 
action on the gas would be stronger for the less massive 
component. In this scenario, it is not surprising that we 
observe a clear offset between the centroid of the X-ray 
residuals (Fig. 18.2(1 and the galaxy density peak around 
BG2, the latter being in advance with respect to the gas 
in the direction of motion. This would be expected if the 
galaxies follow the dark matter and the gas behaves as a 
collisional fluid (see also lF04|> . 

Further evidence for our interpretation comes from the 
"hot spot" centred at (J2000) RA = 22 h 49 m 35.8 s , Dec = 
-64 o 30'25'.'35 (discussed in Sect. 133). We can speculate 
that this is the reverse shock, which can weakly be seen 
in the lower middle panel of Ricker & Sarazin's Fig. 7. 

9. Summary and conclusions 

The main results of our detailed analysis of the XMM- 
Newton observation of the merging cluster A3921, at z — 
0.094 can be summarised as follows: 

— The X-ray morphology of the cluster is elliptical, with 
a pronounced extension toward the NW, in the direc- 
tion of the second and third-brightest galaxies (BG2 
and BG3). The centroid of the X-ray emission is offset 
from the brightest cluster galaxy (BG1) by 17" to the 
NW. 

— A 2D /3-model fit excluding all of the extension to the 
NW shows that the ellipticity of the main cluster is 
e = 0.24. Subtraction of this model from the cluster 
image reveals a large residual structure in the direction 
of the extension. This residual contains both diffuse 
emission and emission from BG2 and BG3. There is 
further significant residual emission at the centre of 
the main cluster. 

— The global temperature estimate (kT = 4.4±0.07 keV) 
is in agreement with previous GINGA and ROSAT re- 
sults. However, the fit is poor (reduced \ 2 = 1-34) and 
there are large residuals around 1 keV, suggesting the 
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Q 25 




Right Ascension 

Fig. 9. X-ray residuals above a 2D-/3-model as in Figure^ 
overlaid on the red sequence galaxy isodensity map with 
a cut of 18 in the I-band magnitude ffrom lF04l The figure 
shows that the secondary X-ray peak is located behind the 
galaxy overdensity in the supposed direction of motion. 
See text for more details. [see the electronic version of the 
paper for a colour figure] 



presence of significant amounts of multi-temperature 
gas. 

The temperature map shows a distinctly nonuniform 
temperature distribution. The temperature structure 
of the Western side is dominated by a "bar" at a sig- 
nificantly higher temperature. The hot bar is aligned 
roughly SE-NW. The Eastern side of the cluster shows 
much less temperature variation. 

A mass profile, calculated under the usual assump- 
tions of hydrostatic equilibrium and spherical symme- 
try from the surface brightness and temperature pro- 
files of the relaxed-looking Eastern side, is not well fit- 
ted by any of the standard profiles from numerical sim- 
ulations. The best-fitting NFW profile yields xt > 4- 
On the other hand, the optical and X-ray mass esti- 
mates agree statistically at a radius of 0.35hj^ Mpc. 
Emission from the brightest BG2 and BG3 is clearly 
visible in both the raw image and in the residual after 
subtraction of the 2D /3-model. The X-ray temperature 
of BG3 is rather high kT = I.I^q "|, and we speculate 
that it may be the centre of a small group of galaxies. 
The X-ray emission from BG2 may come from a weak 
AGN. 



Comparison of the X-ray results with numerical simu- 
lations leads us to conclude that in A3921 a small subclus- 
ter associated with BG2 is undergoing an off-axis merger 



with a larger main cluster. The merger age appears to 
be t < 0.5 Gyr. This is in excellent agreement with the 
optically-derived results of lF04l 

This analysis of the XMM-Newton observation of 
A3921 has allowed us to completely revise the conven- 
tional view of this interesting cluster, which appears cur- 
rently to be the best example of an X-ray observed large- 
scale temperature asymmetry. 

The present observations illustrate the great advan- 
tage to be gained in the study of merging clusters from 
combining observations in several different wavelengths as 
well as with numerical simulations. Further spectroscopic 
observations will give deeper insights into the distribution 
of galaxies in the subcluster. At the same time, the up- 
coming Chandra observation of this cluster should allow 
detailed investigation of the central regions. 
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Appendix A: Analysis of systematic effects in 
temperature determination 

In this Appendix we are interested in assessing possible 
systematic effects on the global properties of the cluster. 

A.l. Nh variations 

In order to have better constraints on temperature when 
building the temperature map, we first extracted a spec- 
trum in the region of the cluster less perturbed by the 
interaction (i.e., by excluding the residuals to the West; 
this is the main region as shown in Fig. |3J). We fitted 
this spectrum with both an absorbed mekal model and 
an absorbed APEC model, where temperature, metallicity 3 
and normalisations (emission measures) were the free pa- 
rameters of the fit. The Galactic column density was also 
left as a free parameter. However, our previous experience 
with XMM-Newton data suggests that it is prudent to 
check for a possible soft excess, which can lead to differ- 
ent temperature values in MOS and pn because of their 
different sensitivities at low energy. We performed several 
tests and concluded that in order for the fitted MOS and 
pn temperatures to agree within the 90 per cent errors the 
channels below 0.3 keV for MOS and below 0.5 keV for 
the pn camera should be excluded. 

The best-fitting results obtained with mekal and 
APEC are statistically equivalent. The N-g. value used for 
the temperature map computation is obtained by simul- 
taneously fitting the MOS spectra only (between 0.3 and 
10.0 keV) with the absorbed mekal model. This gives A^h 
= 2.191q 4 2 x 10 20 cm -2 (90 per cent errors for one sig- 
nificant parameter), which is ~ 25 per cent lower than 
the Galactic value given by Dickey & Lockman (1990), 
but in excellent agreement with the PSPC result of Nh = 
(2.2 ± 0.2) x 10 20 cm" 2 (Arnaud et al. H55j)|> . 

We then investigated possible variations of the Nh in 
the field of view by inspecting the 100/im map at the posi- 
tion of A3921. Adopting the relation between the infrared 
emission at 100/^m and the column density of Boulanger 
et al. (1996), we found values of between 1.93xl0 20 cm~ 2 
and 2.48xl0 20 cm -2 , i.e., a maximum variation of 12 per 
cent. These values are consistent within the errors with 
our best-fitting X-ray value. In view of the good agree- 
ment between all of these measures, throughout the pre- 
sented work, the absorption was fixed to the best-fitting 
value derived from the MOS spectral fit. 

3 We adopted the abundances from Grevesse & Sauval 
(1998). 
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A. 2. Global analysis 

We extracted a spectrum in a large circle of radius 11', 
corresponding to ~ 0.65 r2oo, the global region in Fig. El 
The spectra of the three cameras were then fitted simulta- 
neously with several absorbed thermal models (see Table 
rP|l in the [0.3/0.5 -10.0] keV range (MOS/pn). The free 
parameters of the fits were the temperature and normal- 
isation (emission measure). We also investigated whether 
constraints could be put on abundances of individual ele- 
ments, which also help in evaluating possible systematics 
due to specific element variations. We first fitted with all 
elements as free parameters. However, fitting a vmekal 
or a vapec model to the high S/N global spectrum, we 
were only able to constrain the relative abundances of O, 
Ne, Si, Fe, and Ca. Thus in the fits described below we 
tie all unconstrained elements to Fe. The Galactic column 
density was fixed to 2.19xl0 20 cm~ 2 , as discussed above. 

The best- fitting results are listed in Table fA~T1 Errors 
in the Table are quoted at la for one significant param- 
eter. All model fits give best-fitting temperature, abun- 
dance and normalisation values which agree within their 
90 per cent errors. In most cases (the temperature for ex- 
ample), the best-fitting values agree within the 1 a errors. 
In particular, we note that the derived temperature is not 
sensitive to abundance variations. The best- fitting single- 
temperature model is vapec, although the reduced x 2 = 
1.3 is clearly not satisfactory. 

With the single temperature fit to the global spectrum, 
we observe large residuals corresponding to the Fe L mul- 
tiplet (see the top panel of Fig. [3J, reflecting the exis- 
tence of multi-temperature gas as observed in the tem- 
perature map (Fig. . The largest residuals are centred 
at ~ 0.8keV, and thus we might expect to observe at 
least another component at ~ 0.6 keV (Belsole et al. 2001, 
Bohringer et al. 2002). We thus fitted the spectrum with 
two-temperature thermal models, obtaining significantly 
better fits in all models, as shown in Table lA~TJl . 

In our two-temperature analysis, we observe the fol- 
lowing: 

— we indeed find the 0.6 keV component as expected, 
although its normalisation is about a factor of 10 lower 
than the higher temperature component; 

— the best fit is obtained with two mekal rather than 
two APEC models, which is the opposite of what we 
obtained for the one-temperature analysis above; 

— in both one and two-temperature model fits, the APEC 
model gives higher Fe abundances and lower temper- 
atures than the mekal model (although we note once 
more that all results are statistically in agreement 
within la). This is likely to be a systematic effect (al- 
though a different effect from that found by Buote et 
al. 120031 for example). However, because of the result- 
ing x 2 for the IT and 2T models, it is difficult to esti- 
mate to what extent the temperature/metallicity vari- 
ation is due to fitting a multi-temperature gas with a 
singlc-tcmpcrature model, or due to differences in the 
particular (mekal, APEC) plasma code; 



— the dominant temperature and the Fe and Si abun- 
dances are higher in the two-temperature fits as com- 
pared to the single-temperature fits; the opposite is 
true for O, Ne and Ca. These effects are known as the 
Fe bias and the Si/S bias (see Buote 120001 for details), 
and they underline the complexity of estimating metal- 
licity variations in the presence of a multi-temperature 
gas; 

— in the case of single VMEKAL and vapec model fits, 
the relative abundances of O, Ne, Si and Ca are dif- 
ferent at 90 per cent confidence from the relative Fe 
abundance. Adopting the best-fitting two-temperature 
model, however, the relative abundance of all elements 
are consistent within 1 a with the relative abundance 
of Fe. This result is also partially related to the Fe 
and Si/S biases discussed in Buote (2000). In view of 
this result, and given the asymetric temperature dis- 
tribution of the global region, it is difficult to conclude 
if the relative abundances of the elements above are 
truly different. 

These tests suggest that: 

1. variations in individual abundances have a negligible 
effect on the temperature determination; 

2. neither temperature nor abundance results are 
strongly dependent on the chosen model. 

The existence of the low-temperature model is not re- 
lated to the central galaxy, as discussed in Sect. 14. ll This 
supports our hypothesis concerning the influence of the 
low-temperature component on the relative abundance of 
elements other than Fe. In light of the above we con- 
clude that single element metallicity variation is beyond 
the scope of this paper and anyway is highly uncertain. 
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Table A.l. Best-fitting results of the global spectral analysis. 



single temperature 



Parameter a 



mekal 



apec 

f U.Ob 



vmekal 

QTU5 



vapec 



kT (keV) 

Fe b 

O 

Ne 

Si 

Ca 

Norm c (xl0~ 2 ) 
X 2 /d.o.f. 



4.38(±0.04) 
0.26(±0.02) 



1 86+ - 02 
1551.1/1140 



4.36 



0.30(±0.02) 



1.86(±0.01) 
1549.9/1140 



4.40 
0.26(±0.02) 
0.50(±0.10) 

0.531°;}! 
0.14(±0.06) 

1 i4+°' 31 

1 - i *-0.31 
i OO+0.03 
r.oo_ 02 

1523.2/1136 



4.39(±0.06) 
0.31(±0.03) 



0.771°;}! 
0.82l°;i 4 4 
0.14(±0.07) 

1.79(±0.02) 
1511.9/1136 



two temperatures 



mekal+mekal apec+apec vmekal+mekal vapec+apec 



Parameter 



-FToTTTTT- 

o.ro_ 10 



-0.1(1 

-0.09 

n 44+ - 03 



f0.ll 

J -0.11 
n 30+0 03 

21+ 012 
25+°' 14 

40+°' 08 
0.80lg; 3 3 8 8 

1 71 +0.02 

± -0.02 

fiS + °'° 3 

U - DO -0.03 

16+ 006 

u - lo -0.04 

16+°'° 4 

u - lo -0.04 

1252.6/1132 



5.08lHr 

(-0.03 
-0.03 
j +0.23 
-0.21 
,+0.18 
-0.18 

U - 4O -0.09 
87+0-40 
U - 8 '-0.39 

7n +o.o3 

1 - ' u -0.02 

fi7+ - 04 

u - D '-0.03 

14+009 

-0.03 

1fi+ 004 
1261.9/1132 



kTl (keV) 

Fe b 

O 

Ne 

Si 

Ca 

Norml c (xl0~ 2 ) 
kT2 (keV) 
Fe b 

Norm2 c (xl0~ 2 ) 
xVd-o.f. 



5.09 



,18- 



0.39 



+0.03 



0.44Z 
0.28^ 
0.52;! 



-,+0.02 
J -0.02 
-+0.02 
-0.02 
7 +0.07 
' -0.04 

14+0- 04 
1256.1/1136 



1.70] 
0.65^ 
0.17^ 



0.02 
0.02 
0.03 
0.03 
0.08 
-0.04 

1 cj+0' 04 
1263.7/1136 



1.69] 
0.68^ 
0.17 4 



a Errors are quoted at lcr for one interesting parameter 

b Abundances are the standard photospheric abundances from Grevesse & Sauval 
(1998). For the mekal and APEC models Fe means total abundances. 

c The normalisation is in units of [10~ 14 /(47t(Da x (1 + z)) 2 )] cm~ 5 , where Da 
is the angular size distance to the source (cm), as defined in the mekal model, in 
XSPEC. 



